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Abstract. We have observed a sample of <5 Scuti stars dis¬ 
covered by the ROTSE collaboration in 1999 with Super- 
LOTIS in order to characterize changes in their pulsation 
periods over a time baseline of roughly three years. Achiev¬ 
ing these goals required the creation of an automated 
astrometric and photometric data reduction pipeline for 
the Super-LOTIS camera. Applying this pipeline to data 
from a June 2002 observing campaign, we detect pulsa¬ 
tions in 18 objects, and find that in two cases the periods 
have changed significantly over the three years between 
the ROTSE and Super-LOTIS observations. Since theory 
predicts that evolutionary period changes should be quite 
small, sources of non-evolutionary period changes due to 
the interactions of pulsations modes are discussed. 

Key words: 


1. Introduction 


It is estimated that ~50% of stars located in the lower 
classical instability strip exhibit 6 Scuti pulsations. Most 
pulsate in a number of non-radial p-modes, others in fun¬ 
damental and/or overtone radial modes, and some possi¬ 
bly in g-modes. The periods of these stars are generally 
between 0.014d and 0.25d, with amplitudes up to 1™0. 
There are a number of distinct subgroups of the <5 Scuti 
pulsating stars. In all cases it is the standard k opacity 
mechanism which excites the pulsations, but the stars lo¬ 
cated in this region of the instability strip may be Popu¬ 
lation I, evolved Population II, or massive stars evolving 
through the instability strip. A very thorough description 
of the characteristics of these stars, and their different 


subclasses, is given by Breger (2000). 

The radial pulsation of a 6 Scuti star is related to the 
stellar density through the simple equation of pulsation 
P^fp = Q- Observed changes in the pulsation frequen¬ 
cies of these stars therefore provide important information 
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about stellar structure. Since stellar mass is conserved on 
the time scale of our observations, finding period changes 
provides information about changing stellar radius. De¬ 
tailed models of these changes in stellar structure, and the 
resulting changes in pulsation frequency, have been carried 
out by Breger & Pamyatnykh ( 1998 ) who predicted that 
pulsation frequencies should be stable or slowly increas¬ 
ing as some of these types of stars evolve away from the 
Zero Age Main Sequence (ZAMS). The theoretically ex¬ 
pected exceptions would be the rare pre-MS pulsators, of 
which very few are known, and the Population II stars. 
Breger & Pamyatnykh (1998) point out that in practice 
what is often observed is not a slow evolution of pulsa¬ 
tion frequency, but abrupt changes more drastic th an can 
be explained by the current models. While Szeidl (2000) 
provides examples where period changes have been proven 
false by careful re-examination of data, many examples of 
well established period changes also exist. These observed 
changes are not necessarily even in the predicted sense, 
i.e. some are of increasing frequency, and there are exam¬ 
ples from the literature where the period changes are an 
order of magnitude larger than expected. The predicted, 
and sometimes observed, period changes in J Scuti stars 
are generally given in units of P~ 1 dP/dt « 10~ 7 y _1 , an 
effect on the order of 10~ 8 over the 3 year baseline stud¬ 
ied here. Breger & Pamyatnykh (|l998| ) point out that the 
observed period changes in Population II stars are often 
abrupt and up to an order of magnitude larger than those 
in Population I stars and that stars just below the Main 
Sequence may demonstrate changes up to two orders of 
magnitude larger than Population I stars. 

A small gradient in the period of pulsation requires a 
long baseline of observations of at least a decade to detect 
with certainty. For example, several decades of data for 
the star EH Lib has been analyzed by a number of groups 
(Mahdy & Szeidl (|1980),Yang et al. (1992), and Agerer 


& Huebscher ( 1997[)) who reached different conclusions as 


to subtle changes in the period of pulsation of this star. 
Obviously, systematic effects, and effects related to the 
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method of data analysis, can be problematic when trying 
to identify such a subtle change. Light travel time effects 
due to binarity, and errors with the O-C analysis, or time 
of maxima counting, may lead to such ambiguities in the 
changing pulsation period. Here, large changes in periods 
of pulsation over a comparatively short, three year, time 
scale are investigated in order to compare to theoretical 
predictions. 

2. Sample Stars 

A sample of 18 S Scuti stars was selected from the on¬ 
line database of the Robotic Optical Transient Search 
Experiment (ROTSE) project. The ROTSE team discov¬ 
ered many new variable stars, as describe by Akerlof et 
al. ( |2000|) , and has made the data publicly available^]. 
Stars from the ROTSE database were chosen, as opposed 
to stars from the standard catalogs of S Scuti stars, be¬ 
cause the ROTSE stars were readily observable at the 
time of year when the instrument for this study was avail¬ 
able, and were generally fainter than other samples of S 
Scuti stars, thus better suited to observation by our in¬ 
strument. Several months of well sampled ROTSE obser¬ 
vations could be easily retrieved for each star from the 
online database. Since the ROTSE data were taken with¬ 
out a standard photometric filter, the standard S Scuti 
subclass definitions, based on V band photometric am¬ 
plitudes, are difficult to apply. A selection was made for 
amplitude < O'M and P < 0.25d, and only stars with 
the highest ROTSE light curve quality designation were 
considered. The stars selected from the ROTSE database 
could possibly fall into four categories: (1) Population I p- 
and/or g-mode pulsators, (2) Population I radial funda¬ 
mental and/or overtone pulsators, generally called High 
Amplitude Delta Scutis because they have larger photo¬ 
metric amplitudes, (3) Population II stars, generally called 
SX Phe stars, or (4) pre-MS and evolved massive stars. 
The sample stars and their properties as reported in the 
ROTSE database are given in Table 1. 


day, and is capable of collecting «600 individual 30 sec¬ 
ond images during an eight hour observing run. Since one 
of the primary functions of this instrument is to search for 
the prompt optical emission associated with gamma ray 
bursts, it is able to slew very quickly if a satellite trigger 
is received. As a result, the observations of variable stars 
were occasionally interrupted to chase a burst. 

Images of the S Scuti stars, taken with a Bessel R band 
filter, were typically 20 to 30 seconds in length, with two 
images of the same field taken consecutively. Each star 
was visited typically six times per night, resulting in 10 to 
20 individual data points for each star for each observing 
night. Since the Super-LOTIS CCD camera is thermo- 
electrically cooled, there is a non-negligible dark current 
which must be removed with a dark image. More than 
50 dark images were taken each night, several for each 
unique exposure duration. At several points during each 
night a set of exposures of the random field at Zenith was 
taken for the purpose of constructing a super sky flat by 
median filtering the images of the randomly distributed 
stars. While useful for the study of the periods of 5 Scuti 
stars, these data were also useful for a thorough charac¬ 
terization of the Super-LOTIS instrument. Since the CCD 
camera is located at prime focus of a fast Newtonian opti¬ 
cal system, problems with focus stability and optical aber¬ 
rations, such as coma, were common. Mechanical malfunc¬ 
tions with the drive system and wind induced vibration, 
leading to poorly guided images, as well as electronic prob¬ 
lems with the CCD leading to noisy images, were also 
identified. The 12-bit CCD camera itself has seen many 
years of use, and will be replaced with a device with higher 
quantum efficiency and better read noise characteristics. 
The current CCD has a high gain of 20 e~/ADU and an 
approximate read noise of 70e“. While the data still pro¬ 
duced good results despite these hardware problems, the 
system is constantly being improved. Less data loss, and 
higher quality images, may be expected with the installa¬ 
tion of a newly purchased CCD camera and improvements 
to the mount which will occur during 2003. 


3. Observations 


Observations of the target stars were carried out dur¬ 
ing May and June, 2002, with the Super Livermore Op¬ 
tical Transient Imaging System (Super-LOTIS) robotic 
telescope located at the Steward Observatory site on 
Kitt Peak, Arizon a. T he tel escope system is described by 
Park et al.(1997 and [2002 ). A 2048 2 pixel CCD detec¬ 
tor is located at prime focus of a 60cm f/3.5 Boiler & 
Chivens Newtonian reflector, resulting in a field of view 
of 0.8° 2 and a plate scale of 1.5 ,/ /pixel. Typical seeing at 
the Super-LOTIS site was 2" — 3", so the stellar PSFs 
were generally under-sampled. The system functions au¬ 
tonomously, with a list of targets being submitted each 


1 http://umaxpl.physics.lsa.umich.edu/~mackay/ 
rsvl / rsvl_home.htm 


4. Data Reduction 

A pipeline for the reduction of images from the Super- 
LOTIS telescope was devised for use with this data 
and, hopefully, for use with future projects as well. This 
pipeline was based on codes written in Interactive Data 
Language(IDL), and includes code from public IDL li¬ 
braries, such as the Goddard library]^ The process was 
broken up into three steps: image reduction, astrometry, 
and photometry. First, master dark images were created. 
The dark current was typically measured to be 0.2 ADU/s, 
on top of a bias level of 182 ADU. The original intention 
was to make flat-field images by a median filter of images 
of the random star fields taken throughout the night at 90° 

2 http://idlastro.gsfc.nasa.gov/homepage.html 
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Table 1. The 6 Scuti Star Sample. The listed periods are those determined by the ROTSE team and reported in the 
the ROTSE database. 


Star 1 

RA°(J2000) 

DEC°(J2000) 

Mag 

P(d) 

AMag 

ROTSE 9470 

253.90332 

52.379080 

13.40 

0.1613 

0.39 

ROTSE 4861 

264.90536 

50.200581 

12.45 

0.1761 

0.35 

ROTSE 3956 

194.17738 

23.152025 

13.54 

0.1654 

0.33 

ROTSE 3554 

205.30734 

31.790440 

13.03 

0.1318 

0.31 

ROTSE 1526 

210.27316 

24.704525 

12.24 

0.1939 

0.10 

ROTSE 3124 

214.22528 

23.878872 

13.06 

0.1772 

0.21 

ROTSE 3709 

215.48952 

23.449043 

13.35 

0.1848 

0.14 

ROTSE 7388 

217.66783 

27.224258 

14.23 

0.1744 

0.29 

ROTSE 4837 

221.33072 

35.466972 

13.77 

0.1707 

0.30 

ROTSE 706° 

231.02898 

36.866928 

10.99 

0.1041 

0.33 

ROTSE 9134 

232.84311 

35.881794 

14.32 

0.1784 

0.46 

ROTSE 917 

243.83421 

35.707275 

11.10 

0.1805 

0.28 

ROTSE 2573 

240.62502 

37.560284 

12.31 

0.1794 

0.10 

ROTSE 3232 

242.53889 

35.958557 

12.68 

0.1594 

0.23 

ROTSE 2275 

243.61665 

30.529236 

12.25 

0.1769 

0.10 

ROTSE 2592 

237.44629 

23.915735 

12.45 

0.1110 

0.12 

ROTSE 2926 

242.08601 

28.212255 

12.66 

0.1750 

0.11 

ROTSE 2986 

240.19507 

24.261307 

12.59 

0.1965 

0.13 


(1) These data were collected from the ROTSE variable star databas e, available at 
http://umaxpl. physics, lsa.umich.edu/~mackay/rsvl/rsvl_home.htm 
(a) The known variable star YZ Boo 


altitude. In practice, it was determined that with the expo¬ 
sure time that had been allotted there were not enough sky 
counts to produce a statistically meaningful flat-field. In 
order to investigate the effects of using unflattened data, 
a number of images taken occasionally by Super-LOTIS 
during the morning and evening twilight were carefully 
studied. Since the CCD camera was removed for mainte¬ 
nance on several occasions, it was decided not to attempt 
to apply these few test flats to data taken on other nights. 
For objects within the central 90% of the CCD chip area, 
the contribution to photometric errors from the lack of 
flattening is expected to be < 1%. This is taken into ac¬ 
count during photometric calibration when only stars near 
the center of the CCD chip are utilized for determination 
of the photometric zero point of each image. 

A significant portion of the data, up to 25%, were 
found to be unusable for some reason related to hardware 
or weather conditions. A series of tests were implemented 
in order to reject these poor images at the very onset of 
processing. A requirement was set that there be at least 
100 5cr point source detections in each image. Since data 
were often taken under patchy clouds, this helped to re¬ 
move images taken during very poor transparency or those 
taken too far into the morning twilight. A test was run to 
check the gradient of the sky value across the chip from 
edge to edge, and a cut set at a gradient less than 10%. 
This helped to identify images taken too close to the moon. 
Similarly, a test was run to determine the gradient of the 
sky values in a set of annuli moving from the center of 


the image toward the edge. This test easily identified a 
distinctive problem caused by the occasional shutter mal¬ 
function. Lastly, poor tracking was identified by calculat¬ 
ing the autocorrelation of each image. The characteristic 
’hiccup’ caused by a tracking irregularity in the Right As¬ 
cension gear or mount vibration induced by wind, which 
is easily identifiable by visual inspection, was identified in 
the data pipeline as a strong image autocorrelation well 
outside the normal image seeing disk of « 5 pixels. The 
autocorrelation was also found useful for identifying the 
occasional image with a severely distorted PSF due to 
poor focus. This test was done using FFTs of the images 
and was found to be fast enough as to require only a very 
small portion of the total data pipeline CPU time. 

Once poor quality data were rejected, a World Co¬ 
ordinate System (WCS) astrometric solution was added 
to the FITS header of each remaining hie. Approximate 
center-of-held coordinates are logged by the telescope con¬ 
trol system at the time the exposures are taken, so these 
estimates were used as a starting point for an accurate 
WCS solution. The pointing of the Super-LOTIS system 
was typically better than 6', with an average pointing er¬ 
ror of 5.7' over 1000 pointings. Since the held of view is 
so large, targets were seldom, if ever, not placed on the 
CCD chip, but the pointing accuracy is expected to im¬ 
prove with future upgrades to the telescope mount. The 
astrometric solutions were determined by cross correlat¬ 
ing the positions of point sources within each image with 
the positions of stars in the Guide Star Catalog II. Point 
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sources in the images were identified with an IDL adapta¬ 
tion of the Stetson (1987) DAOPHOT algorithms written 


by W.Landsman and available in the Goddard IDL library. 
There were typically 200-300 matches between image stars 
and catalog stars for each image. The cross correlation of 
the catalog and source positions was maximized allowing 
translation of the image point source positions, a rotation 
of up to 5°, and a low order polynomial term in the radial 
direction. The resulting coefficients were stored in a GSSS 
astrometric structure, of the type used for the HST guide 
star survey images, and added to the image FITS header. 
The typical residuals of the astrometric fits were < 0.5 
pixel. Since the stellar PSFs were under-sampled, the ac¬ 
curacy of the fits was not as high as would be expected 
from the application of t his m ethod to other data sets, 
for example Hogg et al. ( [2001| ), which is followed closely 
here. Parts of the Super-LOTIS pipeline astrometric code 
were drawn from the SP SS P T telescope data pipeline de¬ 
scribed by Hogg et al. ( 2001 ), but here the requirements 
for astrometric precision are much less stringent, « lpixel. 
The failure rate for the astrometric solutions to images 
passing the image quality tests was less than 10%. 

Since the Super-LOTIS images were not crowded with 
stars, the reduced images were photometered using sim¬ 
ple aperture photometry with a radius, matched to the 
average seeing, of 6 pixels and a sky annulus of radii 10 
and 18 pixels. These values were determined empirically 
after testing a variety of radius and annulus combinations. 
Sub-pixel centroiding of the aperture was carried out prior 
to photometering in order to compensate for errors in the 
conversion from Right Ascension and Declination to pixel 
coordinates, which were found to be on the order of 1 pixel. 
Since the stellar PSFs were often very distorted, an accu¬ 
rate centroid was very important. This was done using a 
center-of-mass estimation. In each data image, every star 
amounting to a 5cr detection above the noise in the sky 
background was photometered and its flux, chip position, 
Right Ascension, Declination, and nearby sky value were 
stored, along with the Julian date of the observation, in a 
FITS binary table written to accompany each data image. 
The fluxes were converted to magnitudes by comparison 
to the Guide Star Catalog II (GSCII) on-line database^]. 
The stars photometered in each image were matched to 
stars from the GSCII catalog and their F band photo¬ 
graphic catalog magnitudes recorded. Both F band and 
R band are red bandpasses, and though not identical, are 
very similar. While each individual catalog magnitude has 
a reported error of up to 0™5, a minimum of 100 of these 
values were averaged together in order to substantially re¬ 
duce errors in the photometric zero points. The GSCII 
contains more accurate photometry for bright stars taken 
from other catalogs. For each matched star, a rough zero 
point was calculated with the simple formula 


mj p = —2.5 * log (Flux) + zpt 
3 http://www-gsss.stsci.edu/support/data_access.htm 


( 1 ) 


where mj p is the catalog magnitude from the GSCII, Flux 
is the observed flux, in ADU s _1 , and zpt is the rough pho¬ 
tometric zero point. At least 100 such rough zero points 
were averaged together to create a final zero point for each 
image. Prior to averaging, the rough zero points from the 
brightest and faintest 10% of the stars were removed, and 
the final average was calculated on the remaining values 
with 3(7 clipping. The errors in these catalog-based mag¬ 
nitudes are dominated by the errors in the GSCII magni¬ 
tudes, but were found to be accurate to « 5%. 

Since in many situations accurate relative photome¬ 
try is more important than absolute photometry, the in¬ 
strumental magnitudes were further refined following the 
methods of Everett & Howell ( [2001 ). A small set of stars 
brighter than MTO, and close enough to the target star 
that, based on known pointing accuracy statistics, all of 
these stars would appear in each exposure, were chosen to 
act as the reference for the calculation of differential mag¬ 
nitudes. Assuming that each of the reference stars were 
intrinsically constant in brightness, small corrections were 
made to the instrumental magnitudes following the refer¬ 
ence stars method of Everett & Howell ( 2001 , eqn.(2)). In 
this method the differential magnitudes are calibrated us¬ 
ing an average over N reference stars weighted by the the¬ 
oretical variance of each observation, in magnitudes. This 
estimate was based on the read noise of the CCD, the pho¬ 
ton noise, dark current, as well as noise associated with 
the quantization of the electrons into ADU. The equation 
for er* is taken from Everett & Howell ( [2001 , eqn.(l)), but 
terms for dark current and quantization are added accord¬ 
ing to Howell (2000, eqn.(4.4)) 


1 (1QC7 \/F*~\-Cpi x (F s ky + R 2 + Fdark-\-g 2 *p) 
(J * l.UOO ( Tp 


( 2 ) 


where A* is the aperture flux from the star in e~, g is 
the CCD gain in e~/ADU, F s k y is the average sky back¬ 
ground value in e~ /pixel, Fdark is the average dark cur¬ 
rent in e~ /pixel, R is the CCD read noise in e _ , p is a 
constant relat ed to the quantization noise as explained by 
Howell ( 2000| ), and C P i X relates the number of pixels in 
the aperture to the number of pixels in the annulus via 


Upix — XL a p er \ 1 T (jl a p er /Tl anr /)\ 


( 3 ) 


The accuracy of the differential photometry was com¬ 
pared to the theoretically expected variance in the mag¬ 
nitude estimations derived from the equations for er*. The 
results are shown in Fig.l as function of magnitude, with 
the heavy line presenting the expected level of photometric 
accuracy, due largely to the high read noise of the CCD. 
By plotting the measured variance of a random sample of 
stars, all assumed to be constant, it is shown that the level 
of photometric accuracy achieved is very close to the theo¬ 
retical limit. It was found that for a star of average bright¬ 
ness 12U0 the differential magnitude values produced with 
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this method are accurate to « 2%. A photometric calibra¬ 
tion where a number of stars in a single frame are utilized 
to produce differential magnitudes, sometimes called ’en¬ 
semble photometry’, has been shown to produce highly 
accurate results by several groups, including Everett & 


Howell (2001). This method helps to reduce systematic ef¬ 
fects, as well as effects due to color dependent extinction, 
resulting in accurate differential magnitudes. While the 
best photometric accuracy was achieved when an ensem¬ 
ble of stars was used to produce differential magnitudes, 
the option exists in the pipeline to derive only the less 
accurate results from the GSCII using a larger number of 
stars without requiring that the all the stars be in every 
frame. This is a necessary modification in order to quickly 
produce accurate magnitudes for targets of opportunity, 
such as gamma ray burst optical flashes, when a set of en¬ 
semble stars can not be selected prior to the reductions, a 
large error box must be searched, and photometric results 
are required as quickly as possible. 

5. Period Determination 

With magnitudes for the sample of S Scuti stars produced 
by the photometric pipeline, the periods of variation could 
be determined and compared to the values found from the 
data taken during the 1999 ROTSE observations. Since 
the data are unevenly sampled, a periodogram method 
was selected for analyzing periodicity in the light curves. 
The effectiveness of the periodograms for identifying peri¬ 
odicity in data sampled in a manner similar to the ROTSE 
or Super-LOTIS data has been de mons trated many times, 
i nclud ing Schwarzenberg-Czerny ( 1996 ) and Eyer & Blake 
( 2002 ). Here an implementation of the Lomb-Scargle nor¬ 
malized periodogram, a technique based on the work of 
Lornb ( 197C| ) and Press & Rybicki ( 1989 ), is utilized in 
an IDL code based on a publicly available routine writ¬ 
ten by J. Wilms^.The statistics of various types of pe¬ 
riodograms are treated in great detail by Schwarzenberg- 
Czerny ( |1998|) , in which a common inaccuracy in the stan¬ 
dard implementation of the Lomb periodogram is pointed 
out. In the limiting case where the number of data points 
in the time series is very large, the statistical distribution 
of the periodogram of a white noise vector is well approxi¬ 
mated by a y 2 distribution, but in general the distribution 
is actually a j3 function. In cases were the number of data 
points is less than « 200, the /3 function provides much 
more accurate false alarm probabilities for detected sig¬ 
nals. The full distribution function, and its limiting case, 
is given by 


l-P(z) = (l-^)i 


(4) 


The Lomb method was applied to the publicly avail¬ 
able ROTSE light curves in an attempt to recover the 

4 http://astro.uni-tuebingen.de/software/idl/ 
aitlib/timing/scargle.pro 


periods for the sample stars published in the ROTSE 
database. The period search mesh was set so that an aver¬ 
age resolution in period of 10“' was achieved. Heliocentric 
corrections were made to the times of observations prior to 
analysis with the Lomb-Scargle algorithm. When weighted 
by errors, the periods reported by ROTSE and those de¬ 
termined from analyzing the ROTSE data with the Lomb- 
Scargle algorithm were identical. The published ROTSE 
periods were determined usi ng a cubic B-spline method 
described by Akerlof et al. ( 1994 ). This method is com¬ 
putationally slower than the Lomb algorithm, but is more 
statistically robust for light curves which are not well ap¬ 
proximated by sinusoids. Since the light curves of S Scuti 
stars are basically sinusoidal, the Lomb method is thought 
to be preferable for this application. An example of a peri¬ 
odogram for the star ROTSE 706 is given in Fig.3. While 
the statistical properties of det ected signals have received 
much attention, Akerlof et al. ( 1994 ) points out that fre¬ 
quency error estimation is rarely discussed in the litera¬ 
ture. A simple Monte Carlo method was utilized here to 
estimate the errors in the detected periods. Following Ak¬ 
erlof et al. ( 200C| , Fig.3), the photometric errors reported 
in the ROTSE database were augmented to match the 
published error distribution as a function of magnitude. 
Systematic effects were found to increase the ROTSE pho¬ 
tometric errors from the theoretically derived limit 0.5% 
at 10“0 to ~ 2%. For 200 trials a Gaussian random vec¬ 
tor, of sigma set to match the mean expected photomet¬ 
ric variance of the individual measurements was added 
to the photometric data, processed with the Lomb algo¬ 
rithm, and the resulting highest power frequencies were 
tabulated. The standard deviation of these frequencies was 
taken to be the la error on the originally determined pe¬ 
riod. Periods and accompanying errors for the ROTSE and 
Super-LOTIS photometry using the Lomb algorithm and 
Monte Carlo error estimation are given in Table 2. The sig¬ 
nificance of the period differences are also given, in units 
of the quadrature added errors for the individual period 
estimations. Inspection of our derived period changes, and 
comparison with our Monte Carlo-derived period errors in 
Table 2, may appear to indicate that our period errors are 
overestimated (with reduced chi-squared of <<1 for the 
dataset). However, we note that the Monte Carlo tech¬ 
nique we use will, even in the best circumstances, inflate 
derived errors by a factor of \[2 by adding its “false” er¬ 
rors to the intrinsic errors of the data. Moreover, in our 
numerical investigations of the period-fitting process, we 
found that the period uncertainties are often limited not 
by the shape of the global minimum at the true period, 
but rather by the existence of nearby local minima corre¬ 
sponding to quite distinct period solutions. In this case, 
the traditional Gaussian characterization of the uncertain¬ 
ties is inadequate, and a conservative approach - such as 
we have adopted here - will appear to have overestimated 
the errors in those cases where the true global minima 
have (in fact) been derived. 
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Table 2. Lomb-Scargle periods and Monte Carlo error estimates for the ROTSE data and the Super-LOTIS data 


Star 

ROTSE 

P(d) 

cr(d) 

Super-LOTIS 

P(d) 

a(d) 

b 

oT 

<1 

ROTSE 9470 

0.161347 

l.le-5 

0.161343 

4.8e-5 

0.08 

ROTSE 4861 

0.176131 

6.7e-6 

0.176242 

3.1e-5 

3.50 

ROTSE 3956 

0.1654 

0.0033 

0.1654 

0.0047 

0.00 

ROTSE 3554 

0.1318 

0.0022 

0.1318 

0.0022 

0.00 

ROTSE 1526 

0.193927 

1.0e-5 

0.1939 

0.0077 

0.00 

ROTSE 3124 

0.1772 

0.0027 

0.177173 

9.3e-5 

0.01 

ROTSE 3709 

0.185 

0.010 

0.1850 

0.0090 

0.00 

ROTSE 7388 

0.175 

0.019 

0.175 

0.014 

0.00 

ROTSE 4837 

0.171 

0.011 

0.170756 

6.9e-5 

0.02 

ROTSE 706“ 

0.1040889 

7.0e-7 

0.1040453 

3.9e-6 

11.0 

ROTSE 9134 

0.151 

0.037 

0.151 

0.013 

0.00 

ROTSE 917 

0.1805 

0.0056 

0.180344 

1.5e-5 

0.03 

ROTSE 2573 

0.179 

0.018 

0.179 

0.011 

0.00 

ROTSE 3232 

0.159372 

7.6e-6 

0.159 

0.011 

0.03 

ROTSE 2275 

0.177 

0.016 

0.1768 

0.0053 

0.01 

ROTSE 2592 

0.1199 

0.0025 

0.114 

0.015 

0.40 

ROTSE 2926 

0.175 

0.012 

0.185 

0.059 

0.17 

ROTSE 2986 

0.196518 

1.6e-5 

0.1965 

0.0093 

0.05 


(a) The known variable star YZ Boo 


6. Results 


As expected, for the majority of stars in the <5 Scuti sam¬ 
ple, the principle periods of variation were found, within 
errors, not to have changed between the 1999 ROTSE ob¬ 
servations and the Super-LOTIS observations during May 
and June, 2002. However, two stars, ROTSE 4861 and 
ROTSE 706, were found to exhibit statistically significant 
period variations. In the case of ROTSE 4861 a period 
increase of ss 6.3 ■ 10 -4 was detected at the 3.5er level. 
For ROTSE 706 and period decrease of ^ ~ 1.4 • 10~ 3 
was detected at the ll.Oer level. In Fig.2 the ROTSE and 
Super-LOTIS data sets are combined and phased to the 
ROTSE period. It is clear from the folded light curves of 
ROTSE 706 that there is a phase lag between the two data 
sets which is greater than the error in the phase estima¬ 
tion. At the same time, for the majority of the stars in the 
sample the period errors are large enough that the phase 
uncertainties over the three years separating the observa¬ 
tions becomes large, > .25P, meaning that the phase is 
not expected to be coherent between the two data sets. In 
the case of ROTSE 706, where a period change was iden¬ 
tified and the phase errors due to period inaccuracies were 
negligible, an attempt was made to fit the ephemeris from 
the ROTSE and the Super-LOTIS observations by adding 
a dP/dt term of the magnitude of the determined period 
change. Since these observed period changes are large, and 
therefore not evolutionary, the ephemeris were also fit to 
an abrupt, step-function period change. Neither case was 
found the improve the phase residuals for ROTSE 706. 
Since the period changes could be quite complex in na¬ 


ture, it might be expected that these two simple models 
of period change would fail to adequately fit the data. 


7. Discussion 


The magnitude of these period changes are much larger 
than those expected from stellar evolution, with the mod¬ 
els of Breger & Pamyatnykh ( |199£f ) predicting that Pop¬ 
ulation I 6 Scuti periods in the O.ld to 0.2d range are 
essentially constant, P~ 1 dP/dt < 10 _8 year _1 , over the 
duration of the observations considered here. For pre-MS 
stars the predictions increase by a factor of up to 100, 
still at least an order of magnitude smaller than the pe¬ 
riod changes found here. The non-evolutionary sources 
of per iod c hanges are discussed by Breger & Pamyat¬ 
nykh ( 1998 ). Given the amplitude of the observed period 
changes and the short time span over which these changes 
occurred, it seems unlikely that light travel time effects, 
resulting from binarity, could be responsible. Additionally, 
previous searches for period changes in ROTSE 706 have 
not turned up effects of the magnitude found here, indi¬ 
cating that light travel effects due to a periodic binary 
motion are even less likely. An interesting possibility is 
outlined in the work of Moskalik ( 1985| ) who describes the 
possible effect of the non-linear interaction of various pul- 
sational modes within the star that could produce jumps 
in the pulsation period of a 5 Scuti star of the amplitude 
found here. The theoretical predictions for period changes 
in Population II 6 Scuti stars are less well constrained, 
but it is thought that such changes are typified by large 
jumps in period followed by relatively constant behavior. 
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These jumps are observed, for example in the star CY 
Aqr by Powell et al. (1995|), and are usually of the order 
P~ X AP « 10" 6 . 

Further observations of ROTSE 4861 should be un¬ 
dertaken in order to determine its age and stellar proper¬ 
ties. The star ROTSE 706 is a well studied Population I 
High Amplitude 6 Scuti star also known as YZ Boo. Since 
this star has a high photometric amplitude and relatively 
short period, it has been included in many studies of the 
properties of 6 Scuti stars. The long-term stability of the 
period of pulsation of YZ B oo has been studied exten¬ 
sively by Harndy et al. ( 1986 ) who find a period increase of 
P~ 1 dP/dt = 3 • 10 - 8 year -1 over the last several decades. 
Clearly, our result is incongruous with their finding. Since 
this star is a Population I fundamental radial mode pul- 
sato r, as confirmed by Pena et al. ( 199£ ) and Rodriguez et 
al. ( 1996 ), the abrupt period changes predicted for Popu¬ 
lation II stars are not expected to be present. The most 
promising theoretical explanation for a period change of 
the magnitude observed here in a Population I star is 
non-linear mode interaction. Such interactions have been 
proposed as a possible explanation for the observed pe¬ 


riod changes in XX Pyx reported by Handler et al. (1998 


and 200C). Since 6 Scuti stars may have a large number 
of excited modes of all types, and probably many more 
with amplitudes below current photometric detection lim¬ 
its, the calculation of the specific modes responsible for 
the amplitude and period variations is extremely difficult. 
The mode interactions described by Moskalik (1985) are 
based on only three modes, a radial mode interacting with 
two stable g- or p-modes of low order, and especially con¬ 
sidering that the mode growth factor, 7 , is not precisely 
known for many of the modes identified in XX Pyx, the 
problem of fitting a mode interaction model to the obser¬ 
vations is not feasible. While at one time it was thought 
that monoperiodic High Amplitude 6 Scuti (HADS) did 
not show any evidence for types of pulsation other than 
fundamental radial mode, or sometimes an overtone (see 
Rodriguez ( 1996 )), the mode resonance explanation for 
period changes requires nonradial modes to be excited in 
these stars. The subject of low amplitude nonradial and 
high order radial modes in HADS stars is discusses by Gar- 
rido et al. (1996) who find that diverse types of modes are 
often found in HADS stars when the the photometric data 
is of high enough quality. At the same time, the mode reso¬ 
nance hypothesis requires that the dominant, presumably 
radial, mode be unstable, meaning nonadiabatic. Nonadi- 
abtic simulations of stellar models similar to 6 Scuti stars 
are considered by Stellingwerf (1980) who finds that the 
fundamental radial mode of a stellar model similar to YZ 
Boo is unstable with a mode modulation time scale, 7 , 

of the order years. It is certainly plausible that the con¬ 
ditions, namely multiple modes and an unstable domi¬ 
nant mode, could exist in a 5 Scuti star, such as YZ Boo, 
for period changes of the order A P/P ss 10 -4 to occur 
through resonant mode coupling. There are relatively few 


discussions in the literature of observations of large am¬ 
plitude, non-evolutionary period changes. In addition to 
the case already discussed of XX Pyx, period variations in 
VI162 Ori are presented by Arentoft et al. ( [2001 ), and the 
peculiar behavior the one-time monoperiodic HADS star 
RY Lep is discussed by Laney & Schwendiman (2002). In 
both cases, changes in periodicity, similar in magnitude to 
those presented here, are found. As pointed out by Laney 
& Schwendiman, these abrupt changes in period, or “ape¬ 
riodic cycles” may even be typical of variable stars. Since 
they occur rarely, evidence for such behavior is only now 
becoming more abundant as more data than ever is being 
collected by groups all over the world interested in S Scuti 
stars. 

Comparatively little is known about ROTSE 4861, the 
other star for which we find a period change. The colors of 
the star were measured with The 60” Oscar Meyer Tele¬ 
scope on Mt. Palomar and found to be V — I = 1.37, 
almost a magnitude more red than a typical lower insta¬ 
bility strip star. One interesting note about this star is 
that its position, as reported in the ROTSE database and 
re-calculated by us, is located 5" from the X-Ray source 
1RXS J173936.5+50120 7 from the ROSAT Bright Source 
Catalog of Voges et al. ( 1999| ). This raises the interesting 
possibility that ROTSE 4861 is a new pre-MS, or Herbig 
Ae/Be, star displaying 5 Scuti variability. Only nine such 
stars have been identified, but the discovery, and exten¬ 
sive study, of such objects may prove useful for determin¬ 
ing the time-scales of pre-MS evolution and the structure 
of young stars (see Kurtz (2002)). Further observations of 
ROTSE 4861 are required to determine the nature of this 
star, but Zinnecker & Preibisch (1994) found that about 
half of the Herbig Ae/Be stars they studied were X-Ray 
bright. If found to be a pre-MS star, the large amplitude 
and comparative regularity of ROTSE 4861 may make it 
the best testbed yet found for evolutionary theories. It is 
predicted by Breger & Pamyatnykh ( 1998| ) that for pre- 
MS stars the observed period derivative is negative, which 
it is here, but the amplitude of the change found here is 
about 80 times larger than what is predicted for a pre-MS 
star with a pulsational period of 0.17d. 


8. Conclusions 

Tested with imaging data for a set of 18 S Scuti stars, 
a newly designed data reduction pipeline for the Super- 
LOTIS telescope is found to be efficient and accurate. 
These photometric data were used to compare the cur¬ 
rent periods of these stars to periods determined three 
years ago by the ROTSE project. As expected, the ma¬ 
jority of the stars show the same periodicity in the ROTSE 
data and the Super-LOTIS data. A Monte Carlo simula¬ 
tion was done to estimate the accuracy of the periods de¬ 
termined with a Lomb-Scargle periodogram method, and 
two stars were found to exhibit statistically significant pe¬ 
riod changes. The star ROTSE 706 (=YZ Boo) is a well 
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Fig. 1. The dark line represents the theoretically expected 
photometric accuracy based on equation (2). The squares 
are measurements of the variance in the magnitudes of 
actual stars. 


Fig. 2. The ROTSE and Super-LOTIS data for YZ Boo 
are combined and phased to the same period. The solid 
line with cross points are the ROTSE data and the dashed 
line with box points are the Super-LOTIS data.The sta¬ 
tistically significant phase lag is clearly shown. 


studied monoperiodic high-amplitude 6 Scuti which ap¬ 
pears to have undergone a non-evolutionary period change 
during the time between the two data sets studied here. 
A few other examples of large amplitude period changes 
in 6 Scuti stars were found in the literature. Along with 
these observed period changes was found discussion of a 
possible explanation through the phenomenon of resonant 
mode coupling between a nonadiabatic radial mode and 
two stable nonradial modes. The second period change 
was found in the star ROTSE 4861. Little is known about 
this star, but it is possibly an X-ray source reported in 
the ROSAT BSC. Further observations of this star are 
required to determine if it is a pre-MS or Herbig Ae/Be 
star, in which case it would represent a prime target for 
an astroseismology campaign to learn more about pre-MS 
evolution. 
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Fig. 3. An example of a Lomb-Scargle periodogram gen¬ 
erated from the Super-LOTIS observations of the star YZ 
Boo. 
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